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ABSTRACT 
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Climatological temperature and salinity fields are used to calculate the salinity 
contribution to density and dynamic height fields in the Leeuwin Current System (LCS). 
While the temperature gradient is primarily linear, with warmest water to the north, the 
salinity fields are spatially inhomogenous. A comparison of density fields, calculated 
with constant and variable salinity, shows that, off Western Australia, the density field is 
primarily determined by temperature. Off Southern Australia, the density field is 
dependent on warm and salty (subtropical) and fresh and cold (sub-Antarctic) water 
masses. While the dynamic height fields, calculated with constant and variable salinity, 
show similar flow patterns off Western Australia, different flow patterns are found off 
Southern Australia. 

In addition to the analysis of climatological fields, a primitive equation ocean 
model is used to investigate the role of salinity in the formation of currents and eddies in 
the LCS. Two identical ocean models, one with a climatological salinity field and the 
other with no horizontal salinity gradients, are run and compared with each other. Despite 
the model runs being initialized with similar temperature distributions, there are relatively 
large temperature and density differences in the Southern Australian region, due to the 
advection of water masses by the Leeuwin Current. 

Based on the climatological analyses and the results of the model experiments, it 
is concluded that, descriptively and dynamically, both temperature and salinity are 
essential to accurately characterize the large-scale circulation of the LCS. 
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I. INTRODUCTION 



At a given pressure, the density of sea water depends on temperature and salinity. 
Except for high-latitude regions and near sources and sinks of salt, temperature is usually 
considered to vary more than salinity in the upper ocean and is generally considered to be 
the major source of density variations (Pickard and Emery, 1990). 

Recent dynamic height analyses by Batteen et al. (1995) have shown that the 
distribution of salinity in the California Current System (CCS) can be important in 
defining the large-scale circulation of the CCS. In particular, the mean variability of 
salinity has been shown to be responsible for a significant equatorward component along 
the coast of California and a strong offshore component adjacent to Baja. Batteen et al. 
(1995) concluded that, descriptively and dynamically, both temperature and salinity are 
essential to accurately characterize the large-scale structure of the CCS. 

Besides the CCS, other eastern boundary current (EBC) regions, such as the 
Leeuwin Current System (LCS), could show significant changes if salinity as well as 
temperature effects are considered. The Leeuwin Current is an anomalous EBC which 
flows poleward (against the prevailing wind direction) along the Western Australian 
coast, down to Cape Leeuwin (See Figure 1 for geographical locations), and swings 
eastward extending as far east as the Great Australian Bight (e.g., Cresswell and Golding, 
1980). There is general agreement (e.g., Godfrey and Ridgway, 1985) that the Leeuwin 
Current is generated by a meridional pressure gradient, which overwhelms the opposing 
wind stress. The source for the Leeuwin Current water is predominantly geostrophic 
inflow from the west (e.g., McCreary et al., 1986; Thompson, 1987) and is augmented by 
a source from the North West Shelf (e.g., Gentilli, 1972), possibly having its origin in the 
Pacific Ocean (Godfrey and Ridgway, 1985). 

Since the Leeuwin Current is a density driven current, both temperature and 
salinity could contribute to its development. While the temperature gradient is primarily 
linear, with warmest water to the north, the salinity distribution is spatially 
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inhomogeneous (Figure 2). Off of Western Australia low salinity water is present to the 
north and northwest, high salinity to the southwest, and low salinity to the south. 
Depending on the time of year, these water masses could contribute to the development 
of the Leeuwin Current, subsurface currents, and eddies. 

Here we ask the question whether salinity variability is an important consideration 
for defining the large-scale circulation of the LCS. To quantify the importance of salinity 
in the LCS, we calculate its contribution to density and to dynamic height in the LCS 
(Chapter II). After incorporating salinity into an ocean model (Chapter III, Section A), 
two identical ocean models, one with a climatological salinity field and the other with no 
temporal or spatial salinity gradients, are run and compared with each other (Chapter III 
Section B). Although the results with variable salinity included show similar features to 
the results with constant salinity off Western Australia, there are important quantitative 
differences in the Cape Leeuwin and Southern Australian regions. The low and high 
salinity tags in the LCS are found to be particularly useful for identifying the 
contributions of the different water masses to the development of the Leeuwin Current 
and eddies. 
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II. ANALYSIS OF HYDROGRAPHIC AND DYNAMIC HEIGHT FIELDS 



A. ANALYSIS OF ANNUAL HYDROGRAPHIC AND DYNAMIC HEIGHT 
FIELDS 

Two major water masses occur for all or part of the year off Western Australia 
(Figure 2): a warm (tropical) low salinity water mass from the North West Shelf (e.g.. 
Weaver and Middleton, 1989) possibly having its origins in the Pacific Ocean (Godfrey and 
Ridgway, 1985), and a cooler (subtropical), higher saline water mass called the 'Western 
Australian Current' which is influenced by the presence of a large semi-permanent trough 
extending in a northeastward direction (Andrews, 1 977). 

The seas just south of Cape Leeuwin and eastward to the entrance (near Esperance 
in Figure 2) to the Great Australian Bight axe influenced by several water types: 1) 
subtropical water from west of Western Australia carried in a branch of the Western 
Australian Current in summer (Andrews, 1977), 2) tropical water from north of Western 
Australia brought by the Leeuwin Current in autumn and winter (Rochford, 1969; Cresswell 
and Golding, 1980), and 3) local Sub-Antarctic water characterized by its coldness and low 
salinity (Rochford, 1986). 

To explore the water masses of this area, the climatological data sets of Levitus et 
al. (1994) and Levitus and Boyer (1994) were employed. The data are gridded with one- 
degree horizontal resolution at standard depth levels. These data sets, which include 
annual and monthly mean temperature and salinity fields, were used to compute annual 
and seasonal density and dynamic height fields. 

Before presenting the annual and monthly climatological temperature, salinity, 
density, and dynamic height fields, we describe how the calculations for density and 
dynamic height were obtained. Following Batteen et al. (1995), the density is computed 
using the following approximate equation of state for seawater: 

P = Pot 1 " «(T -r 0 ) + p(S - S o )], 
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where a and (3 are the expansion and contraction coefficients for temperature T and 
salinity S, respectively, p is density, S 0 is the reference salinity, and T 0 is the reference 
temperature. This equation can estimate the effect of salinity for a given region if typical 
values of a, P, T 0 , and S 0 are used. For the LCS region, a is chosen to be 2. 5X10 -4 (°K)'', 
P is 7.5XKT 4 , T 0 is 278.2 °K and S 0 is 35.5 psu (practical salinity units). 

To assess the effect of salinity on density in the LCS, we calculate density in the 
above equation for variable and for constant salinity. Calculating density for variable 
salinity means that the density in the above equation is computed at a given level with 
climatological values of T and S. Calculating density for constant salinity means that the 
density in the above equation has S replaced by S 0 , so that the above equation reduces to 

p = p 0 [\-a(T-T 0 )]. 

To assess the effect of salinity gradients on the large-scale ocean circulation of the 
LCS, dynamic height fields are calculated with and without the inclusion of salinity 
variations. The dynamic height anomaly AO is calculated by integrating the specific 
volume anomaly 5 over depth and comparing it with A0 35 5TP (where the subscripts 35.5, 
T, and P stand for salinity 35.5 psu, temperature, and pressure, respectively), calculated 
by integrating 5 X . Comparisons between horizontal gradients of AO and A0 35 5TP are used 
to provide insight concerning the importance of salinity variations to the flow field. The 
reference level used in this study is 1000 m depth. 

The climatological annual mean temperature fields are shown at 10, 100, and 200 
m depth in Plate la, and at 400 m depth in Figure 3a. At all depths shown, the 
temperature distribution is primarily zonal, with the isotherms running roughly east-west. 
Maximum temperatures are in the north and decrease to the south. The temperature 
differences from north to south are ~ 1 1°C, 9°C, 7°C, and 2°C at 10, 100, 200, and 400 m 
depth, respectively. 

The climatological annual mean salinity fields are shown at 10, 100, and 200 m 
depth in Plate lb, and at 400 m depth in Figure 3b. The common features of the 10 and 
100 m depth plots can be clearly seen: low salinity (tropical) water to the north and 
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northwest, high salinity (subtropical) water to the southwest, and low salinity (sub- 
Antarctic) water to the south. The maximum salinity for 10 m depth is centered at ~ 32°S, 
while the maximum salinity for 100 m depth is centered at ~ 28°S. The salinity gradient 
weakens from 10 to 100 m depth. The 35.5 salinity contour extends all along Southern 
Australia at 10 m depth, but only as far east as Albany at 100 m depth. At 200 m depth, 
the salinity pattern is similar to the zonal temperature pattern with maximum salinity in 
the north and minimum salinity in the south. At this depth, the 35.5 salinity contour 
extends only as far south as Fremantle (~ 32°S). At 400 m depth, the salinity difference is 
only 0.4 psu from north to south. 

The annual density fields calculated with constant (variable) salinity at 10, 100, 
and 200 m depth are shown in Plate lc (Plate Id). The annual density fields calculated 
with constant (variable) salinity at 400 m depth are shown in Figure 3c (Figure 3d). A 
comparison of the upper levels (i.e., 10 and 100 m depth) in Plate lc and Plate Id shows 
the following. Although, off Western Australia the density is primarily determined by the 
temperature field, the density fields in Plate Id are generally less dense than those in Plate 
lc due to the presence of the low salinity (tropical) water mass. South of ~ 34°S, the 
effect of salinity on density becomes significant. In particular, off Southern Australia, the 
density is influenced by both the warm and salty (subtropical) water mass and the fresh 
and cold (sub-Antarctic) water mass. At 200 m depth, the density in Plate Id is generally 
greater than that in Plate lc due to the presence of high salinity throughout the region. At 
400 m depth, the density field in both figures is relatively uniform throughout the region. 

Figures 4 and 5 show the annual dynamic height patterns at 10, 100, 200, and 400 
m depth for constant salinity, i.e., A 0 3 5 5 T p, and for variable salinity, i.e., Ad>, 
respectively. In both figures, the upper level (i.e., 10 and 100 m depth) dynamic height 
contours tend to be perpendicular to the west coast of Australia, which is consistent with 
onshore geostrophic flow toward Western Australia. Between 34°S and 36°S, while the 
upper level dynamic height fields in Figure 4 tend to be zonally distributed, those in 
Figure 5 tend to vary in both the meridional and zonal directions, due to the influence of 
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the salty (subtropical) water mass offshore. The lower level (i.e., 200 and 400 m depth) 
dynamic height patterns in both figures show equatorward flow off Western Australia, 
which is consistent with the depths at which the equatorward undercurrent has been 
observed (e.g., Thompson, 1984). 

B. ANALYSIS OF SEASONAL HYDROGRAPHIC AND DYNAMIC HEIGHT 
FIELDS 

Since the LCS has a strong seasonal (e.g., the Leeuwin Current is strongest from 
March to August) as well as annual cycle, the monthly mean temperature and salinity 
climatological fields are examined to evaluate the effects of the seasonal evolution of 
salinity on the density and dynamic height fields. The annual mean temperature and 
salinity fields showed that the water mass characteristics of the LCS are very similar at 10 
and 100 m depth. Since the seasonal cycle has essentially the same vertical pattern in the 
upper layers, we present only the 10 m depth seasonal fields here. 

The same three water masses described earlier are discussed here to help identify 
how salinity seasonally influences the density and dynamic height fields. These water 
masses are: the warm and fresh tropical water mass off northwestern Australia, the warm 
and salty subtropical water mass off southwestern Australia, and the cold and fresh water 
mass off Southern Australia. 

Plate 2a shows the monthly mean temperature fields in February, May, August, 
and November. As the November and February plots show, the temperature gradient off 
Western Australia during the austral spring and summer has a 45 degree slope from 
northeast to southwest, which provides the source for an onshore geostrophic current (i.e., 
the Leeuwin Current). As the May and August plots show, the slope becomes even 
steeper during the austral fall and winter, which is consistent with the Leeuwin Current 
being strongest during this period. Off Southern Australia, except for the advection of the 
Leeuwin Current in the vicinity of Cape Leeuwin during the austral fall and winter, the 
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temperature pattern follows a typical seasonal cycle of warmest (coldest) temperature off 
the coast in summer (winter). 

Plate 2b shows the monthly mean salinity fields in February, May, August, and 
November. In February, the salinity field shows a zonal distribution of fresh water in the 
north, which becomes saltier poleward to - 32°S. There is an offshore region of 
maximum salinity of - 35.8 from - 28°S to 35°S, centered at - 32°S. South of - 35°S, 
fresher water is found. In May, the salinity gradient has a 45 degree slope off Western 
Australia. The maximum salinity has shifted southward to - 34°S. In August, the fresh 
water off northwest Australia has moved farther southward, while the maximum salinity 
has moved farther offshore. In November, off northwest Australia, the fresh water has 
been replaced by more saline water. The maximum salinity has retreated equatorward to 
- 31°S. 

The density fields shown in Figures 6a and 6b were calculated with constant and 
variable salinity, respectively. In general, throughout the year, the density fields off 
Western and Southern Australia in Figure 6b are less dense than those in Figure 6a, due 
to the presence of low salinity (tropical and sub-Antarctic) water masses. In the 
southwestern Australia region, the density fields in Figure 6b are more dense than those 
in Figure 6a because of the presence of the salty (sub-tropical) water mass. 

Figures 7a and 7b show the seasonal dynamic height patterns for constant and 
variable salinity, respectively. In both figures, throughout the year, off Western Australia, 
the dynamic height flow pattern is consistent with onshore geostrophic flow toward 
Western Australia. In contrast, south of 36°S, while the dynamic height fields in Figure 
7a tend to be zonally distributed throughout the year, those in Figure 7b tend to be vary in 
both the meridional and zonal directions, due to the influence of the salty (subtropical) 
water mass offshore. 

To quantify the seasonal effect of salinity on the dynamic height field, the 
dynamic height field at several offshore and nearshore locations were calculated. Figure 8 
shows the meridionally averaged dynamic height field at 108.5°E from 1) 22°S to 29°S, 
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2) 30°S to 35°S, and 3) 36°S to 40°S. In region 1, the dynamic height with variable 
salinity (solid line) is similar to that with constant salinity (dashed line). In the other two 
regions, the dynamic height fields are dependent on both temperature and salinity from ~ 
May to August in region 2 and for most of the year in region 3. In region 2, the seasonal 
dependence on both temperature and salinity coincides with the time frame that the 
Leeuwin Current is strongest. In region 3, the seasonal salinity effect is consistent with 
the time that the subtropical salty water moves poleward to ~ 40°S. 

The dynamic height fields calculated at several nearshore locations are shown in 
Figure 9. In the nearshore regions 1 and 2, which correspond to coastal locations off 
Western Australia, because the salinity gradient is not as strong as in the offshore regions, 
the salinity effect is not very significant in these regions. In contrast, in the nearshore 
region 3, which corresponds to a coastal location off Southern Australia, the dynamic 
height is dependent on both temperature and salinity throughout most of the year, 
indicating that in this region salinity plays a very important role. 
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III. EFFECT OF SALINITY ON AN OCEAN MODEL OF THE LCS 



A. THE MODEL 

1. Model Description 

The numerical model in this study was originally used for a coarse resolution, 
closed basin by Haney (1974), and later adapted for eddy-resolving, limited EBC regions 
with open borders on the northern, western, eastern, and southern boundaries by Batteen 
et al. (1989, 1996) and Batteen (1996). The limited area EBC model is multi-level, uses 
non-adiabatic primitive equations on a beta-plane, and has both baroclinic and barotropic 
velocity components. The model is based on the hydrostatic, Boussinesq, and rigid lid 
approximations. The governing equations are as follows: 



- 1 dp d 2 u 


(1) 


-1 cb , d 2 v 

i v + K M -^ 

p 0 dy J M M dz 2 


(2) 


du dv dw 
H + — 0 

dx. dy dz 


(3) 


ii 

i 


(4) 


p=/3 0 [i-a(r-r 0 )+/?(S-s 0 )] 


(5) 


dt H H dz 1 


(6) 


dS , d 2 S 


(7) 



In the above equations, t is time, ( x,y,z ) is a right-handed Cartesian coordinate 
system with x pointing toward shore, y alongshore, and z upward. The corresponding 
velocity components are (u,v,w), T is temperature, S is salinity, p is density, and p is 
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pressure. Table 1 provides a list of other symbols found in the model equations, as well as 
values of constants used throughout the study. 

For the finite differencing, a space-staggered B-scheme (Arakawa and Lamb, 
1977) is used in the horizontal. Batteen and Han (1981) have shown that this scheme is 
appropriate when the grid spacing is approximately on the same order as, or less than, the 
Rossby radius of deformation, which meets the criteria of this study. The horizontal grid 
spacing is 14 km in the north-south direction and 1 1 km in the east- west direction, while 
the internal Rossby radius of deformation is ~ 30 km. In the vertical, the 10 layers are 
separated by constant z-levels of 13, 46, 98, 182, 316, 529, 870, 1416, 2283, and 3656 m. 
This spacing scheme concentrates more on the upper, dynamically active part of the 
ocean, above the thermocline. 

The model domain (Figure 2) is a rectangular region encompassing the west and 
south coasts of Australia, from ~ 22.5°S to 40°S (1792 km alongshore), and from 
~107.5°E to 122.5°E (1408 km cross-shore). The coastal boundaries of the model domain 
are closed, and have both the tangential and normal components of velocity set to zero. 
Bottom topography has been omitted to focus on the roles played by thermal and saline 
forcing. The constant depth used in the model is 4500 m. A modified version of the 
radiation boundary conditions of Camerlengo and O'Brien (1980) is used for the open 
ocean domain boundaries to the north, south, east and west. Some spatial smoothing is 
applied in the vicinity of the open boundaries. 

The model uses biharmonic lateral heat and momentum diffusion with the same 
choice of coefficients (i.e., 2.0 x 10 I7 cm 4 s"’) as in Batteen et al. (1989). Holland (1978) 
showed that the highly scale-selective biharmonic diffusion acts predominantly on 
submesoscales, while Holland and Batteen (1986) found that baroclinic mesoscale 
processes can be damped by Laplacian lateral heat diffusion. As a result, the use of 
biharmonic lateral diffusion should allow mesoscale eddy generation via barotropic 
(horizontal shear) and/or baroclinic (vertical shear) instability mechanisms. As in Batteen 
et al. (1989), weak (0.5 cm 2 s’ 1 ) vertical eddy viscosities and conductivities are used. 
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Bottom stress is parameterized by a simplified quadratic drag law (Weatherly, 1972), as 
in Batteen et al. (1989). 

2. Method of Solution 

Equations (1) through (7) comprise a closed system of seven scalar equations and 
seven unknowns, u, v, w, p, p , T, and S. The variables, u, v, T, and S are prognostic 
variables whose time rates of change are predicted from (1), (3), (6) and (7), respectively. 
Although the diagnostic variables w, p, and p can be determined from (3), (4), and (5), 
respectively, there are additional constraints imposed on p and w by the choice of the 
rigid lid boundary conditions. The vertically integrated pressure can no longer be 
obtained by integrating the hydrostatic equation (4) for the free surface. Further, the 
vertically-integrated horizontal velocity is constrained to be non-divergent, i.e.. 




( 8 ) 



which is obtained by integrating (3) and applying the vertical boundary conditions. 

For any quantity q, let its vertical average be denoted by q and its departure 

(vertical shear) by q' . From (8) the vertical mean flow can then be described by a 
streamfunction if / , such that: 



u = 



1 dy/ 



1 dy/ 

HHc 



( 9 ) 

( 10 ) 



The streamfunction y/ is predicted from the vorticity equation, which is derived 
from the vertical average of (1) and (2). Applying the curl operator to the vertical average 
(1) and (2), and using (9) and (10), the vorticity equation is 
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where G and F represent the collected contributions of the nonlinear and viscous terms in 
equations (1) and (2). 

The vorticity equation (11) is solved by obtaining an updated value of C, by 
application of the leapfrog (or every 1 1 time steps, the Euler-backward) time-differencing 
scheme. The associated value of if/ can then be obtained from: 
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( 12 ) 



{dy 2 ) dx dx dy dy 

which is an elliptic equation. A solution to (12) is fully prescribed by specifying the 
values of \f/ on the open and closed boundaries of the model domain. Currently, to solve 
(12), the model uses an elliptic solver when there are no variations in coastline geometry 
and/or topography, and successive over-relaxation techniques when there are variations in 
coastline geometry and/or topography. 

The vertical shear current (u\ v') is predicted from (1) and (2) after subtracting the 
vertical mean. The results are: 



a 2 u' 



du' - 1 dp' , . 

— = — + fv'-A M V*u’ + K M , 

dt p 0 ck J M M a 2 



— T 

+ F-F- 



PoH ’ 



dv' -\ dp' , d 2 v' — r" 



(13) 



(14) 



dt p 0 dy J M M cz l p 0 H 

In (13) and (14), p', which represents the departure of the pressure from the 
vertical average, is , using (4), expressed in terms of p as: 
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( 15 ) 



P' = [pgde ~^i H {pg ds )dz 

where £ is a dummy variable representing the vertical coordinate. 

The method of solution consists of predicting y. /, u' ,v' ,T, and S from (1 1), 

(12), (13), (14), (6) and (7), respectively. The total current is then obtained by adding the 
vertical shear part to the vertical average part, after the latter is obtained from (// using 

(9) and (10). The diagnostics p, w, and p' are obtained from (5), (8), and (15), 
respectively. 

3. Forcing Conditions and Experimental Design 

Previous experiments by Batteen and Rutherford (1990) and Batteen et al. (1992) 
investigated the role of thermal forcing in generating the Leeuwin Current and eddies off 
Western Australia with a model domain which was a closed eastern boundary with open 
borders to the north, south, and west. The model was initialized with climatological 
thermal forcing from Levitus (1982), and the ocean was then allowed to geostrophically 
adjust in the absence of external forcing. 

In a more recent study by Batteen et al. (1996), thermal forcing conditions from 
Levitus et al. (1994) and Levitus and Boyer (1994) were used. In addition, the forcing 
was continuously applied at the western boundary, and the model domain was expanded 
to include the Southern Australian region, an area rich with observations of currents and 
eddies, but lacking in eddy-resolving modeling studies of the region. Model 
developments included the opening up of the eastern boundary in the vicinity of 
Esperance in Southern Australia, and the incorporation of an irregular coastline into the 
model. 

In this study, to explore the additional effect of salinity on the LCS, both annual 
and seasonal salinity forcing conditions from Levitus et al. (1994) and Levitus and Boyer 
(1994) have been incorporated into the model. The annual (seasonal) temperature and 
salinity forcing conditions for the upper five levels (top level), which are assumed to be 
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zonally homogenous, are shown in Figure 10 (Figure 11). For both annual and seasonal 
forcing conditions, since the lower five levels do not exhibit much horizontal variation, 
they are assumed to be constant for each level. The temperature values used for levels 6 
to 10 are 9.53°C, 6.03°C, 3.24°C, 2.19°C, and 1.27°C. The salinity constant used for the 
lower five levels is 34.7. 

The solid, dashed, and dotted lines plotted in Figure 1 1 correspond to latitudinal 
regions where the three water masses discussed earlier are found. Regardless of the 
region, the temperature conditions (Figure 11a) indicate that there is, as expected, a 
temperature maximum in March (austral summer) and a temperature minimum in 
September (austral winter) throughout the whole region. The salinity conditions (Figure 
1 lb) show relatively strong aperiodic (non-seasonal) variations throughout the year for all 
of the regions presented. 

The design of the model experiments (see Table 2) is as follows. Experiments 1 
and 2 study the model response to annual forcing. Experiment 1 uses the variable thermal 
and salinity forcing conditions shown in Figures 10a and 10b. Experiment 2 uses the 
same variable thermal forcing conditions but has salinity set to a constant such that S is 
replaced by S 0 in equation (5). Experiments 3 and 4 study the model response to seasonal 
forcing. Experiment 3 uses the variable thermal and salinity forcing conditions shown in 
Figures 1 la and 1 lb. Experiment 4 uses the same variable thermal forcing conditions but 
has salinity set to a constant such that S is replaced by S 0 in equation (5). 

B. RESULTS FROM MODEL SIMULATIONS 

In each of the experiments, the model integrations start from a state of rest. Once 
a day, the model is updated at the western boundary with annual or seasonal forcing 
conditions. To ensure that the time period when the Leeuwin Current is strongest (~ April 
to September), is incorporated in the model, each model integration is run for a full year. 

Although we are interested in the difference between model runs with variable 
salinity and with constant salinity, we first show that the model produces reasonable 
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results. The main phenomena in the LCS is the density-driven Leeuwin Current, which 
flows poleward off Western Australia and eastward off Southern Australia; an 
equatorward undercurrent found below the Leeuwin Current off Western Australia; and 
eddies. Figures 12 and 13 show that the Leeuwin Current with velocities of ~ 30-100 
cm/s, an equatorward undercurrent with core velocities of ~ 5-15 cm/s, and anticyclonic 
eddies (one just south of Shark Bay, one between Dongara and Fremantle, and one in the 
vicinity of Cape Leeuwin) are reproduced in the model. The results from model runs are 
qualitatively similar and compare well with available observations of the LCS (e.g., 
Boland et al., 1988; Cresswell and Golding, 1980). 

1. Experiments 1 and 2 - Annual Thermal and Salinity Forcing 

Experiments 

The annual forcing model integration with variable salinity (Experiment 1) and 
that with constant salinity (Experiment 2) begin with identical temperature distributions. 
When the two integrations are compared with each other, significant differences are 
found, as is illustrated here for the time period that the Leeuwin Current is strongest ( ~ 
April to September, corresponding to model days 90 to 273). Differences of up to 4°C are 
noted during this time period (Figure 14a), due to the advection of warmer (21°C or 
greater) water in Experiment 1 by the Leeuwin Current between Cliffy Head and Albany 
(compare Figure 15a with 15b). Farther east, between Albany and Esperance, the 
temperature difference of ~ 3°C in Experiment 1 is due to the trapping of warmer water 
nearshore by the Leeuwin Current. 

The time-averaged salinity field (Figure 16) shows that the three major water 
masses discussed earlier (see Figure 2) are reproduced well in the LCS model, i.e., fresh 
(<35.5) water is found off the northern part of Western Australia, salty (>35.5) water 
offshore of Western Australia and Cape Leeuwin, and fresh (<35.5) water to the south. 

The time-averaged density fields (Figures 17a and 17b), along with the density 
difference between Experiments 1 and 2 (Figure 14b) show the following. The density 
differences north of 34°S are relatively small. The water for Experiment 1 off the 
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northern part of Western Australia is less dense than that in Experiment 2, due to the 
combination of warm and fresh water in the region, while offshore of Western Australia, 
the water for Experiment 1 is more dense, due to the combination of salty and warm 
water. 

The largest density differences (Figure 14b) are found between Cape Leeuwin and 
Esperance. The water off Cape Leeuwin in Experiment 1 is more dense than that in 
Experiment 2 due to the saltier water. The water between Cliffy Head and Albany is less 
dense in Experiment 1 than in Experiment 2, due to the much warmer (up to 4°C) water 
in the region. The water east of Albany is more dense in Experiment 1 than in 
Experiment 2, due to the saltier water nearshore which has been advected into the region 
by the Leeuwin Current. 

2. Experiments 3 and 4 - Seasonal Thermal and Salinity Forcing 

Experiments 

Since the LCS has significant temporal variability in salinity (see Figure lib), 
Experiments 3 and 4 are designed to study the model response to seasonal forcing, both 
using the same thermal forcing conditions but with variable (Experiment 3) and constant 
(Experiment 4) salinity. As in the previous experiments, here we focus on the difference 
between the two model runs to identify the salinity effect on density. During the time 
period that the Leeuwin Current is strongest, differences of up to 2°C are noted (Figure 
1 8a) due to the displacement of warm, anticyclonic eddies by the Leeuwin Current in the 
Southern Australian region. 

North of 34°S, the density differences (Figure 18b) show that the water for 
Experiment 3, compared to Experiment 4, is less dense off Western Australia and more 
dense offshore of Western Australia. The water off Cape Leeuwin in Experiment 3 is 
more dense than that in Experiment 4 due to the combination of salty and warm water in 
the region. The water between Cliffy Head and Albany is less dense, due to the 
displacement of warm (up to 2°C), anticyclonic eddies by the Leeuwin Current in the 
region. 
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IV. SUMMARY 



This study used climatological temperature and salinity fields to calculate the 
salinity contribution to density and dynamic height fields in the Leeuwin Current System 
(LCS). While the temperature gradient is primarily linear, with warmest water to the 
north, the salinity field is spatially inhomogeneous with low salinity (tropical) water to 
the north and northwest, high salinity (subtropical) water to the southwest, and low 
salinity (sub-Antarctic) water to the south. A comparison of annual density fields, 
calculated with constant and variable salinity, showed that, off Western Australia, the 
density field was primarily determined by temperature. South of 34°S, the effect of 
salinity on density became significant. In particular, off Southern Australia, the density 
was shown to be influenced by both the warm and salty (subtropical) water mass and the 
fresh and cold (sub-Antarctic) water mass. The seasonal density field patterns were found 
to be similar to the annual density field patterns. 

The annual dynamic height fields, calculated with constant and variable salinity, 
showed similar flow patterns off Western Australia but different flow patterns between 
34°S and 36°S. Throughout the year, the seasonal dynamic height fields calculated with 
constant and variable salinity showed similar flow patterns off Western Australia and 
different flow patterns south of 36°S. 

In addition to the analysis of the climatological temperature and salinity fields, a 
high resolution, multi-level, primitive equation model was used to investigate the role of 
salinity in the formation of currents and eddies in the LCS. After incorporating salinity 
into the model, two identical ocean models, one with a climatological salinity field and 
the other with no horizontal salinity gradients, were run and compared with each other. 

Experiments 1 and 2 used annual thermal and salinity conditions to initialize and 
force the model. Despite the model runs being initialized with the same temperature 
distributions, differences of up to 4°C were noted during the time period that the Leeuwin 
Current was strongest in the Southern Australian region, due to the advection of warmer 
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water by the Leeuwin Current between Cliffy Head and Albany. Large density 
differences were also found between Cape Leeuwin and Esperance, where the salinity 
effect on density became significant. 

Experiments 3 and 4 used seasonal thermal and salinity conditions to initialize 
and force the model. During the time period that the Leeuwin Current was strongest, 
differences of up to 2°C were noted in the Southern Australian region, due to the 
displacement of warm, anticyclonic eddies by the Leeuwin Current nearshore. Large 
density differences were also found between Cape Leeuwin and Albany, where the 
salinity effect on density became significant. 

The results from both the climatological analysis and the model experiments 
strongly suggest that there are important quantitative differences between constant and 
variable salinity calculations, particularly in the Cape Leeuwin and Southern Australian 
regions. The low and high salinity tags in the LCS are found to be particularly useful for 
identifying the contributions of the different water masses to the development of the 
Leeuwin Current and eddies. 
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Figure 1. Box shows the domain of the model for the Leeuwin Current System 
(LCS) off Western and Southern Australia. Domain is bounded by 22.5°S to 40°S, 
107.5°E to 122.5°E. The irregular coastline (solid line) used in the model includes 
Shark Bay (SB), Dongara (D), Fremantle (F), Cape Leeuwin (CL), Cliffy Head 
(CH), Albany (A), and Esperance (E). 
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INDONESIA 




Figure 2. Schematic chart of the mean large-scale circulation in the eastern Indian 
Ocean. The tropical (salinity < 34.5) and subtropical (> 35.8) water masses are 
shaded. The edge of the continental shelf (200 m isobath) is dashed. The Leeuwin 
Current is the southward flow of warm low-salinity water down the Western 
Australian coast, (from Pearce and Cresswell, 1985) 
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Figure 3. Annual temperature (a) and salinity 7 (b) fields at 400 m depth. Annual 
density 7 field with constant (c) and variable (d) salinity 7 at 400 m depth. The contour 
interval is 0.2°C in (a), 0.05 psu in (b) and 0.05 in (c) and (d). 
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Figure 4. Annual dynamic height field with constant salinity at (a) 10, (b) 100, (c) 
200, and (b) 400 m depth. The contour interval is 2 dyn. cm. 
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Figure 5. Annual dynamic height field with variable salinity at (a) 10, (b) 100, (c) 
200, and (d) 400 m depth. The contour interval is 2 dyn. cm. 
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Figure 6. Seasonal density (a t ) fields with constant (a) and variable (b) salinity at 10 
m depth in February, May, August, and November. The contour interval is 0.5 
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Figure 7. Seasonal dynamic height fields with constant (a) and variable (b) salinity 
at 10 m depth in February', May, August, and November. The contour interval is 2 
dyn. cm. 
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Figure 8. Seasonal dynamic height fields meridionally averaged at 108.5°E from (1) 
22.5°S to 29°S; (2) 30°S to 35°S; and (3) 36°S to 40°S. The dashed (solid) line refers 
to dynamic height calculated with constant (variable) salinity . 
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Figure 9. Seasonal dynamic height fields in the nearshore region at (1) 113.5°E, 
which is latitudinally averaged from 22°S to 29°S; (2) 115.5°E, which is latitudinally 
averaged from 30°S to 35°S; (3) and 36.5°S, which is longitudinally averaged from 
116°E to 125°E. The dashed (solid) line refers to dynamic height calculated with 
constant (variable) salinity. 
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Figure 10. The annual temperature (a) and salinity (b) field forcing conditions for 
the upper five levels. Since the lower five levels do not have much latitudinal 
variation, they are assumed to be constant for each level. 



31 




depth=182m 



depth =31 6m 





b 



32 



25.0 - 




Figure 11 Time series plot at 10 m depth of monthly temperature (a) and salinity (b) 
fields used as seasonal forcing in Experiments 3 and 4. The solid lines are for 40°S, 
the dashed lines for 31.5°S, and the dotted lines for 22.5°S. 
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Figure 12. Experiment 1 surface density (a,) field and velocity vectors at day 90. In 
this and the following figures, to avoid clutter, velocity vectors are plotted at ever)' 
fourth grid point in the alongshore and cross-shore directions. The contour interval 
is 0.5. 
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Figure 13. Cross-shore section of velocity field for Experiment 1 at day 90. Solid 
lines indicate equatorward flow, while dashed lines indicate poleward flow. The 
contour interval is 5 cm/s. 
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Figure 14. Difference of time-averaged temperature (a) and density (b) from April to 
September (days 90-273) between annual forcing experiments with variable 
(Experiment 1) and constant (Experiment 2) salinity. The contour interval is 1°C in 
(a) and 0.1 in (b). 
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Figure 15. Time-averaged temperature fields from April to September (days 90-273) 
for (a) Experiment 1 and (b) Experiment 2. The contour interval is 1°C. 
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Figure 16. Experiment 1 time-averaged salinity fields from April to September (days 
90-273). The contour interval is 0.1 psu. 




41 



DEPTH : 13m 
T (DAY) : 90 to 273 



O 

3 



24.0°S 



26.0°5 



28.0°5 



30.0°S 



32.0°S 



34.0°S 



36.0°S 




38.0°S 

1 09.0°E 



25JD' 



1 1 1 .0°E 1 1 3.0°E 

-> 100 . 



115.0°E 

LONGITUDE 



1 

117.0°E 



« k o 
f 

1 19.0°E 




121. 0°E 



Figure 17. Time-averaged density (a t ) fields and velocity vectors from April to 
September (days 90-273) for (a) Experiment 1 and (b) Experiment 2. The contour 
interval is 0.5. 
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Figure 18. Difference of time-averaged temperature (a) and density (a t ) (b) fields 
from April to September (days 90-273) between annual forcing experiments with 
variable (Experiment 3) and constant (Experiment 4) salinity. The contour interval 
is 1°C in (a) and 0.1 in (b). 
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Plate 1. Climatological annual mean temperature (a) and salinity (b) fields at 10, 
100, and 200 m depth. Note that the coastline differentiates with depth. Annual 
density fields with constant (c) variable (d) salinity at 10, 100, and 200 m depth. 
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Plate 2. Climatological seasonal temperature (a) and salinity (b) fields at 10 m 
depth in February, May, August, and November. 
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Table 1. Values of Constants Used in the Model 



Constant 


Value 


Definition 


To 


278.2°K 


Constant Reference Temperature 


So 


34.7 


Constant Reference Salinity 


Po 


1.0276 gm cm 3 


Density of Sea Water At T 0 and S 0 


a 


2.4 x 10* 4 (°K)'‘ 


Thermal Expansion Coefficient 


P 


7.5 x KT 4 


Saline Expansion Coefficient 


K 


10 


Number of Levels In Vertical 


Ax 


1.1 x 10° cm 


Cross-Shore Grid Spacing 


Ay 


1 .4 x 1 0° cm 


Alongshore Grid Spacing 


H 


4.5 x 10 3 cm 


Total Ocean Depth 


At 


800 s 


Time Step 


fo 


0.75 x 10 4 s' 1 


Mean Coriolis Parameter 


g 


980 cm s 2 


Acceleration of Gravity 


A M 


2 x 10 ,; cm 4 s' 1 


Biharmonic Momentum Diffusion Coefficient 


A h 


2 x 10 w cm -4 s’ 1 


Biharmonic Heat Diffusion Coefficient 


Km 


0.5 cm 2 s’ 1 


Vertical Eddy Viscosity 


K h 


0.5 cnf s" 1 


Vertical Eddy Conductivity 
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Table 2. The Design of The Experiments 



Experiment # 


Experiment 


Response To Annual Forcing 


1 


Variable Thermal Forcing with Constant Salinity Conditions 


2 


Variable Thermal and Salinity Forcing Conditions 


Response To Seasonal Forcing 


3 


Variable Thermal Forcing With Constant Salinity Conditions 


4 


Variable Thermal and Salinity Forcing Conditions 
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